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TECHNICAL FIELD 
relates to digital filtering, and in particular to digital 



The present invention 
filter design- 



BACKGROUND 



some digital signal processing algorithms are based on filter design m ti^e 
frequency domain. Typical appUcations are noise suppression algorithms for 
speech enhancement (see [1, 2]) and non-linear processors for echo cancel- 
lation. These filters are often quite long, which means that it may be desir- 
able to perform the convolution (filtering) in the frequency domain, smce this 
is less complex than convolution in the time domain. Since the input signal 
to be filtered is typically much longer than the filter, the desired bnear 
convolution has to be implemented by circular sectioned or block convolu- 
tion (see (31) in order to avoid unacceptable delays and/or complexity. 

A problem with filters designed in the frequency domain is that they are real- 
valued, which leads to a time domain representation in which the peak of the 
filter is spUt between the beginning and end of the filter (this is equivalent to 
a filter that is symmeixic around lag 0. i.e. an non-causal filter). This makes 
the filter unsuitable for circular block convolution, since such a filter will 
generate temporal aliasing. This problem may be mitigated (but not ehmi- 
n^ted) by weighting the data with a window that spans more than one block 
(see [21). However, this will introduce a delay of 1 block, which is undesir- 
able. The split peak also makes the filter unsuitable for time domain convo- 
lution, since the important parts of the filter are at the beginning and end of 
tlie filter. This makes it difficult to approximate the filter with a shorter filter 
t4 reduce the complexity of the time domain convolution. 
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Refer^tei^if 'i^scrft^es a method that starts with an initial reduced length 
discrete frequency filter response (using cither the Barlett or Welch method). 
The reduced length filter response is transformed to the time domain using a 
reduced length IDFT, circularly shifted to obtain a linear phase filter, zero- 
5 padded to obtain an esctended length to avoid temporal aliasing, and trans- 

formed back to the firequcncy domain using an extended length DFT. A 
drawback of this method is that the resolution in frequency is reduced due 
to the initial reduced length filter. 

10 SUMMARY 

An object of the present invention is to provide a filter design method that 
avoids both temporal aliasing and the delay generated by multi-block 
windowing with improved finequency resolution. 

15 

This object is achieved in accordance with the attached claims. 

Briefly, the present invention circularly shifts the filter in the time domain to 
recenter the filter peak and applies a short window around the peak to 
20 extract the essential part of the filter and automatically obtain zero-padding. 

This method has several advantages: 

1. The resulting zero-padded reduced length linear phase filter may be 
used to implement circular block convolution without temporal 

.25 aliasing. 

2. If the window is sufficiently short, the convolution may also be per- 
formed in the time domain. 



30 



3: 



If desired, the obtained linear phase filter may be transformed into a 
minimum phase filter to reduce the algorithmic processing delay. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention, together with further objects and advantages thereof, may best 
be understood by making reference to the following description taken together 
with the accompanyir^ drawings, in which: 

Pig. 1 is a time diagram illustrating a sample sequence x[n]; 

F^. 2 is a time diagram illustrating a filter impulse response h[nl; 

Fig. 3 is a time diagram illustrating a time shifted filter impulse re- 
sponse h[n-l]; 

Fig. 4 is a time diagram illustrating a time shifted filter impulse re- 
sponse h[n-L]; 

Fig. 5 is a time diagram illustrating the convolution of sample secjuence 
x[nl and filter impulse response hin]; 

Fig. 6 is a time diagram illustrating a zero-padded sample sequence 

xfn]; 

Fig. 7 is a time diagram illustrating a zero-padded filter impulse re- 
sponse h[n]; 

Fig. 8 is a time diagram of a longer sample sequence x(n]; 

Fig. 9 is a time diagram of a first zero-padded block of sample sequence 
x{n] in fig. 8 padded with zeroes; 

Fig. 10 is a time diagram of a second zero-p>added block of sample 
sequence x[n] in fig. 8; 

Fig. 11 is a time diagram of a third zero-padded block of sample se- 
quence x[n] in fig. 8; 

Fig. 12 is a time diagram of a diagram of the convolution between said 
first block in fig. 9 and the zero-padded filter impulse response in fig. 7; 

Fig. 13 is a time diagram of a diagram of the convolution between said 
second block in fig. 10 and the zero-padded filter impulse response in fig. 7; 

Fig. 14 is a time diagram of a diagram of the convolution between said 
third block in fig. 11 and the zero-padded filter impulse response in fig. 7; 

Fig. 15 is a frequency diagram illustrating a real-valued filter transfer 
function H[k]; 
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^T^lsTLe diagram of the filter impulse response hN that is 
obteined after transforming the transfer function in fig. 15 to the time domam; 

Fig. 17 is a time diagram of the filter impulse response of fig. 16 after a 
circular shift of N/2 samples; 
5 Fig. 18 is a time diagram illustrating a time window w[nl; 

Fig. 19 is a time diagram of the filter impulse response in fig. 17 after 
multiplication by the window in fig. 18; 

Fig. 20 is a time diagram of the filter impulse response of fig. 19 after a 

circular shift to remove leading zeroes; 
-Lo Fig. 21 is a frequency diagram of the magnitude of the filter transfer 

function Ftkl that is obtained after transforming the filter impulse response of 
fig. 20 to the frequency domain; 

Fig. 22 is a frequency diagram of the phase of the filter transfer function 
F[k] that is obtained after transforming the filter impulse response of fig. 20 to 

15 the frequency domain; 

Fig. 23 is a time diagram of the filter impulse response of fig. 20 after 

transformation to a minimum phase filter; 

Fig. 24 is a fiow chart of an exemplary embodiment of the convolution 
method in accordance with the present invention; and 
2 0 Fig. 25 is a block diagram of an exemplaiy embodiment of a convolution 

apparatus in accordance with the present invention. 

DETAILED DESCRIPTION 

i 25 Since the concepts linear and circular convolution are essential for the 

ptesent invention, these concepts wQl be described in more detail with 
reference to fig. 1-14. 



Fig. 1 is a time diagram iUustrating a short sample sequence x[nl of length L. 
30 This sequence is to be convolved with a filter having an impulse response hln] 

: illustrated in fig. 2. The definition of linear convolution in the time-domain is 

. (assuming a linear time-invariant sjrstem): 
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y[n]^xin}S>h[n]^ Z x[ik]A[n-*) = in this case = 



This expression may be interpreted in the following way: 

5 L Form a set of shifted sequences {h[n], h[n-l], ... , h[n.(L-l)l} as iUus- 

trated in Gg- 2-4, 

2. Multiply each sequence by a corresponding sample xjO], x[lj, x[L-l] 
of sequence x[n]. 

10 

3. Add the scaled sequences to form the convolution as illustrated in fig. 5. 

From the above description and from fig. 5 it is clear that the resulting 
convolution y[nl will have a length of L+M-1 samples. Thus, the filtered 
15 sequence yin] wiU be longer than the original sequence x[nl, which has only L 

samples. 

In the above description it was assumed that the convolution was performed 
in the time domain- However, it is also possible to perform convolution in the 
20 frequency domain by transforming each of the sequences x[nl, h[n] to the 

frequency domain using the Discrete Fourier Transform, DFT (typically 
implemented by the Fast Fourier Transform. FFT), multiply the corresponding 
rr. transforms with each other and perform an Inverse Discrete Fourier Trans- 

form, IDFT (typically implemented by the Inverse Fast Fourier Transform, 
*55 IFFT). back to the time domain to obtain ylnj. Thus, in tiie frequency domain 

. / . the convolution may be expressed as: 



0 30 This method is called circular convolution. In fact this method is often pre- 

: : : ferred, since it reduces the complexity of the convolution, especially for longer 
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filters. However, there are some problems with this approach that have to be 
solved: 

1. Typically the input sequence x[nj and the filter sequence h[n] ore of 
5 different lengths^ but the transfomied sequences have to be of the same 

length to allow multiplication of the transforms. This problem may eas- 
ily be solved by zero-padding the shorter sequence up to the length of 
the longer sequence before the DFT operations. 

10 2. As noted above, the filtered sequence y[n] has a length of L+M-1 sam- 

ples, while the input sequence xln] has a length of only L samples. The 
first M-1 samples of this output sequence are erroneous due to so called 
temporal aliasing caused by a wrap aroiand of the last M-1 samples in 
the output sequence (the last M-1 samples in fig 5 would be added to 

15 the first M-1 samples). Hie result is referred to as circular convolution 

with aliasing. However, what we want is a circular convolution that is 
equivalent to a linear convolution containing all L+M-1 samples of y[n]. 
The proper way to handle this problem is to zero-pad both x[n] and h[nl 
up to at least the length N=L+M-1 of y(n] before performing the DFT, as 

20 illustrated in fig. 6 and 7. 



Fig. 8 is a time diagram of a longer sample sequence x[n]. The longer sequence 
may be filtered in the time donuiin using the definition of linear convolution in 
the time domain given above. Since the filter only influences a limited (M) 
25 input samples at a time, it is possible to calculate the output sample of y[n] 

"on the as the filter h[nl is shifted through the input sequence xln]. Thus, 
only the latest M samples of x[n] are required to calculate the current value of 
yW- However, if the filter is long the method is quite complex. In such cases a 
fi:equency domain method would be preferable. 

In principle it would be possible to filter lor^ sequences in the frequency 
domain in the way described above. However, practical considerations often 
make this approach iinacceptable. For example, consider a sequence of 
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sp^ sa-np... in a »tephony appBcation. speech i, "^^^^ 

^ptog rate of 800O ,amples/s. Thus. 1 second of speech wdl contsan 8000 
Zp^s the con.p.«d.y of «.e neceses^ DFT and rDVT «ould ^ 
r for such ^ se^ces. *no«>er pn>h.em is «»t such an ^pp.^ 

wouW ^ to unacceptable detays. since the entire input se^ence haa ^ be 
collected before the DFT can be calcula«d. A better approach « to u,e block 

convolution in accordance with either the o,erlap-«.ve or the o«rla^^d 
^thod. The overlap-add method win now be described wth reference tn fig. 
8-14. 

m accordance with the overlai^add method the long input sequence xW of fig. 
8 is divided into blocks of length L as iUustxated by fig. 9-11. Each block is 
«ro-padded up to the require! length L^M-l (when. M is the filter length). 
Each zero-padded block is convolved with tHe zero-padded filter (se fig. 7^ m 
the frequency domain as described with reference to fig 6-7.. Tlus 
results in the convolved blocks illustrated in fig. 12-14. Finally, the 
convolved sequences are added to form the output sequence y[n\. It is 
noted that in the overlap regions indicated in fig. 12-14 the final result 
is obtained by adding partial results from two block convolutions. This is 
the reason for the name overlap-add. 

The overlap-save method is similar to the overlap-add method. However, 
in this case the blocks of the input sequence xln] are overlappmg (the 
last samples of one block are identical to the first samples of the next 
block). After frequency domain convolution, the erroneous first M-1 
samples of each filtered block (the samples that contain temporal 
aliasing) are discarded before the blocks are reassembled. 

in some appUcations the filter is determined in the frequency domain. For 
example, in telephony applications noise suppression based on spectral 
subtraction is often used (see [1. 2]). In this case the filter is determined as a 
function H(ca) of frequency: 
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if(a>) = 




where a, 6 are constants and and <i>x(^) are estimates of the 

5 power spectral density of the pure noise and noisy speech, respectively. This 

expression is obtained from the model: 

10 where ^r[n] is the noise signal, x[nl is the noi^ speech signal and ytn] is the 

desired speech signal. An estimate of the desired signal ytn], in which the 
noise has been suppressed, is obtained by applying the filter represented by 
H(u) to the noisy signal x[n]. 

15 Another example of an application in which the filter is determined in the 

frequency domain is a frequency selective non-linear procsessor for echo 
cancellation. In this case the filter is defmed by the function: 



/f(fl>)=:/(l>,(d>).6^(<»)) 

where represent an estimate of the power spectral density of a signal 

xlnl contaminated by echo and ^g(<o) represents an estimate of the power 
spectral density of the echo signal e[n]. The filter is be^d on the model: 



20 



where yln] is the desired speech signal. An estimate of the desired signal y[n], 
in which the echo has been cancelled, is obtained by applying the filter 
represented by H(ca) to the echo contaminated signal xln]. 



14:36 FAX +46 18 153050 AROS P/iCENT AB - PATENTVERKET 3013 

+46 18 153050 



WctPatert-ochTBg-verket 

1999 -09- 0 7 9 

Huvudfoxen Kossan 

in the above examples the filter is described by a real-valued continuous- 
Liuency transfer Action HN- This function is sampled to obtam a d^screte- 
frequency transfer function HlkJ. This step is typical^ performed when «ie 
estimates ar. based on parametric estimation methods. However, it is also 
possible to obtain the discrete-ftequenpy transfer function H[kl direct^, for 
example by using periodogram based estimation methods. An advantage of 
parametric estimation methods is that the estimates typically have lower 
variance than estimates from periodogram based methods. 

Fig 15 is a ft^iuency diagram iUustrating an example of a real-valued discrete- 
frequency filter ttansfer function H[k]. Fig. 16 is a time diagram of the filter 
impulse response h[n] that is obtained after transforming the transfer function 
in fig. 15 to the time domain using an IDFT. As may be seen from fig. 15 this 

filter has some unattractive features: 

1 The filter has the full length N. which means that the alrea^ty calculated 
transfer fimction H[k] of fig. 15 can not be used for frequency domam 
convolution without causing temporal aliasing as described above. 

2. If the convohition is performed in the time domain, the filter is not 
suitable, since its peak is spUt between the beginning and end of the 
filter h[n]. This will introduce a long algorithmic processing delay. 

in accordance with the present invention both these disadvantages may be 
avoided. This will now be explained witii reference to fig. 17-22. 

The first step is to circularly shift the filter h\n] by N/2 samples. Mathemati- 
cally this may be expressed as: 



h[n\ -> h[in -^N/l) mod/i] 
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Fig. 17 is a time diagram of the filter impulse response of fig. 16 after a 
circular shift of N/2 samples (N is assiomed to be an even integer in this 
^cample). 

5 The next step is to apply a window to the shifted filter. Fig. 18 is a time 

diagram illustrating an example of a time window w(n]. In this case the 
window is a Kaiser window (see [5]) having a length of M samples (in the 
example M is assumed to be an odd integer less tiian N). Fig. 19 is a time 
diagram of the filter impulse response in fig. 17 after multiplication hy the 
10 window in fig. 18. The result is a filter u[n] with leading and trailiivg zeroes. 

Preferably the leading zeroes are removed by another circialar shift of D 
samples, where D is the number of leading zeroes. Fig. 20 is a time diagram of 
the filter impulse response of fig. 19 after a circular shift to remove leading 
15 zeroes. The resulting filter fin] will have its first non-zero tap at n-0 and will 

have N-M trailing zeroes. 

The result of the steps performed in fig. 17-20 is illustrated in fig. 21-22. Pig. 
21 is a fi^quency diagram of the magnitude of the filter transfer function Flk] 
20 that is obtained after transforming the filter impulse response fin] of fig. 20 to 

the fi«quency domain. Fig. 22 is a corresponding frequency diagram of the 
phase of the filter transfer fimction F(k]. As may be seen from these f^ures the 
real-valued filter Hlk] has been transformed into a linear phase filter F[k] 
(actually linear phase is obtained already after the first circular shift in fig. 
: '25 17). 



30 



It is noted that the transformed filter fin] in fig. 20 has an effective length of 
only M taps, and that the remaining N-M taps are zero. This fact may be 
exploited in two ways: 

1. The zeit>es may be used to implement linear convolution in the fre- 
quency domain by using circular convolution of blocks of length L in 



■99 09/07 14:37 FAX +46 18 153050 AROS P/ATENT AB PATENTVERKET BlOlS 

-t-4e 18 153050 

Inlet Patent- ochreg-verttd 
1999 -09- 0 7 

11 

Huvudfoxen Kossan 

combination with the overlap-arid method (fig. 8-14) or the overlap-save 
method. The block length L is obtained as: 



10 



15 



20 



2. 



The convolution may be performed in the time domain if the filter length 
M is selected sufficiently small. 



The GSM system will now be considered as an example of this procedure. This 
system already has a natural block length for block convolution, since speech 
is divided into frames of 160 samples. Thus, a natural choice is a convolution 
block length L of 160 samples. This leads to a natural FFT length N of 256 
samples. This results in a maximum filter length M of 97 samples (M<-N- 
L+1). 

For a convolution block length of 256 samples the filter length typically Ues in 
the interval 20-97 filter taps. For a Kaiser window the parameter p typically 
Ues in the interval 2-5. Experiments have shown that a filter length M of 55 
samples obtained fix>m a Kaiser window having a length of 55 samples and a 
parameter p»3 is a good choice, which even allows implementation of the 
convolution in the time domain using a signal processor. 



Fig. 20 illustrates that the linear phase filter f[n] is symmetric around the peak 
at (M-l)/2. This will lead to an algorithmic delay P of half the filter length. This 

='.^ 25 algorithmk: delay may be reduced by transforming the filter f[n] to a minimum 

phase filter. A minimum phase filter may be obtained flnom the linear phase 
filter in fig. 20 by using the method described in 16]. Briefly, this method starts 

*•'•* with the expresmon: 



..30 
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where m are the taps of a linear phase filter that is re-centered (linearty 
shifted) so that it is synunetric around tap 0 instead of tap (M-l)/2 as in fig. 
20. According to the spectral factorization lemma f(z) may also be written as: 

where C is a constant and g(z) is a polynomial 



10 



15 



20 



having all zeroes inside the lanit circle. If g(z) is determined, a minimiim 
phase fUter with the same magnitude response in the frequency domain as 
the linear phase fUter may be obtained as: 

Cgiz)g{z) 

Several methods exist for determining the coefficients of g(z) from the 
coefficients of f(z). For example, g(z) may be determined by sohring the non- 
linear system of equations: 



ffo] 
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25 



using the Newton-Rapson algorithm. Applying the spectral factorization 
lemma to the linear phase fUter in fig. 20 results in a minimum phase filter 
with the same magnitude response (fig. 21) in the frequency domain. In the 
time domain the algorithmic delay P in fig. 22 has been reduced, as illus- 
trated in fig. 23. 
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In the above description of fig. 15-23 a certain processing order has been 
assumed. For example, the circular shift in fig. 17 was performed before the 
windowing in fig- 19. However, this order may be reversed if the window is 
circularly shifted accordingly. Furthermore, the circular shift in fig. 20 may 
5 be eliminated by including it in the circular shift in fig. 17 and circular^ 

shifting the window accordingly. 

In this description a Kaiser window was used. However, all the common 
window types, such as Hamming, Hanning. Barlett windows etc. are 
10 feasible. Even a rectangular window may be used. 

Pig. 24 is a flow chart of an exemplary embodiment of the convolution method 
in accordance with the present invention. This embodiment relates to noise 
suppression and uses fi«quency domain convolution based on the overlap-add 
15 method and minimiam phase filtering. Step SI collects the next input signal 

block. Using this block step S2 determines the transfer function H[k] fi:x>m the 
above described noise suppression algorithm. Step S3 transforms H\k\ to the 
time domain using the IFFT. Step S4 circularly shifts the time domain repre- 
sentation of the filter N/2 samples to produce a linear phase filter. Step S5 
20 applies a window to the linear phase filter to reduce its length. Step S6 

performs another circular shift to remove leading zeroes created by the 
window. Step S7 transforms the resulting reduced length linear phase filter 
into a minimum phase filter. Step S8 transforms the already zero-padded filter 
and the zero-padded input signal block to the firequency domain using the 
■ '25 FFT. Step S9 multipUes the filter transform with the signal block transform to 

• perform the convolution. Step SIO transfonns the result back to the time 

domain using the IFFT. Step Sll combines the current convolved block with 
the previous block in accordance with fig. 12-14. Thereafter the algorithm 
returns to step SI to collect the next input signal block and repeat the 
30 procedure. 



If the convolution is performed in the time domain, steps S8- SI 1 are replaced 
by a time domain convolution step. 
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Fig 25 i« a block diagram of an exemplary embodiment of a convolution 
apparatus in accordance with the present invention. Hiis embodiment is 
suitable for performing the method described in fig. 24. An input signal x[nl is 
forwarded to a buffer 10. which collects and stores a signal block. Usmg the 
input signal block a filter design block calculates the noise suppression filter 
HIk]. An IFFT block 14 transforms the filter to the time domain. A circular 
shift block 26 circularly shifts the time representation of the filter by N/2 
samples. A window block 18 appUes a window to the shifted filter. Another 
circular shift block 20 removes leading zeroes from the reduced length filter. A 
minimum phase btock 22 transforms the filter into a minimum phase filter. 
An FFT block transforms the filter back to the frequency domain. The trans- 
form is forwarded to a multipUcation block 26. The input signal block from 
buffer 10 is also forwarded to a zero-pad block 20. which pads zeroes to the 
block up to the length of L+M-1. The zero-padded btock is transformed to the 
frequency domain by an FFT btock 30. The transform is fonvardcd to the other 
input of multiplication block 26. The convolved signal block from multipUca- 
tion block 30 is transformed back to the time domain in anlFFT block 32. The 
time representation of the convolved block is forwarded to an overlap buffer 34 
and to a combination block 36. Overlap buffer 34 attracts the overlap part 
(the last part) of the received signal block and outputs the stored overlap part 
of the previous convoWed signal block. Combination block 36 adds the overlap 
part of the previous block to the beginning of the current block and outputs 
the remaining part of the current block unchanged. The result is the con- 
volved signal y[n]. 

Typically the blocks in fig. 25 are implemented by one or several microproc- 
essors or micro/ signal processor combinations. 

If the convolution is performed in the time domain, blocks 24- 36 are replaced 
by a time domain convolution block. 
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Rather than repeating the above description, it is noted that a digital fUter 
design apparatus in accordance with an exemplaiy embodiment of the 
present invention wUl include blocks 10-22 in fig. 25. 

5 It will be understood by those skiUed in the art that various modifications 

and changes may be made to the present invention without departure from 
the scope thereof, which is defined by the appended claims. 
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1 . A method of designing a digital filter, including the steps of 
deterniining a real-valued discretc-fixqucncy representation of a desired 

5 distal filter, 

transforming said discrete-frequency representation into a correspond- 
ing discrete-time representation; characteriacd by 

circularly shifting said discrete-time representation; and 

applying a window to said discrete-time representation to produce a 

1 0 zero-padded reduced length filter. 

2. The method of claim 1, characteitaed by circularly shifting said 
reduced length filter to remove leading zeroes. 
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3. The method of claim 1 or 2, ohazaeteriaed in that said real-valued 
discrete-fi^ency representation is formed by a noise suppressing spectral 
subtraction algorithm. 

4. The method of claim 1 or 2, chaiaetexiMd in that said real-valued 
2 0 discrete-frequency representation is fonned by a fiiequency selective non-linear 

algorithm for echo cancellation. 

5. The method of any of the preceding claims, chaiacteriaed in that said 
window is a Kaiser window. 
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6. The method of any of the preceding claims, ehaiacterised by trans- 
forming said reduced length filter into a minimum phase filter. 

7. A digital convolution method, including the steps of 

determining a real-valued discrete-frequency representation of a desired 

digital filter, 

transforming said discrete-frequency representation into a correspond- 
ing discrete-time representation; characterised by 
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circularly shifting said discrete-time representatiDn; 
applying a window to said discrete-time representation to produce a 
zero-padded reduced length filter, and 

convolving an input signal with said zero-padded reduced length filter. 

8. The method of claim 7. cha«cteri«ad by circularly shifting said 
reduced length filter to remove leadix^ zeroes. 

9. The method of claims 7, chaiacteftoed transfonning said reduced 
length filter into a minimum phase filter. 

10. The method of claim 7. 8 or 9. duwacteriaed by performing the convo- 
lution step in the time domain using tiic discrete-time representation of said 
reduced length filter. 

11. The method of claim 7, 8 or 9. cha««teri«ed by performing the convo- 
lution step in the frequency domain by using tiie overlap-add metiiod. 

12. The metiiod of claim 7. 8 or 9, cbaMcteriaed by performing the convo- 
lution step in the frequency domain by using the overlap-save metiaod. 

13. A digital filter design apparatus, including 

means (12) for determining a real-vahied discrete-firequency representa- 
tion of a desired distal filter, 

means (14) for transforming said discrete-frequency representation into 
a corresponding discrete-time representation; characteitaed by 

means (16) for circularly shifting said discrete-time representation; and 

means (18) for applying a window to said discrete-time representation to 
produce a zero-padded reduced length filter. 



14. The apparatus of claim 13. characterised by means (20) for circularly 
shifting said reduced length filter to remove leading zeroes. 
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15. The apparatus of claim 13 or 14, characterised in that said window 
applying means (18) implements a Kaiser window. 

16. The apparatus of any of the preceding claims 13-15, characterisBcd by 
5 means (22) for transforming said reduced length filter into a minimum phase 

filter. 

17. A digital convolution apparatus, including 

means (12) for detemiining a real-valued discrete-frequency representa- 
10 tion of a desired digital filter, 

means (14) for transforming said discrete-frequency representation into 
a corresponding discrete-time representation; chaiacterlzcd by 

means (16) for circularly shifting said discrete-time representation; 
means (18) for applying a window to said discrete-time representation to 
15 produce a zero-padded reduced length filter; and 

means (24-36) for convolving an input signal with said zero-padded re- 
duced length filter. 

18. The apparatus of claim 17, characterised by means (20) for circularly 
2 0 shifdng said reduced length filter to remove leading zeroes. 

19. The apparatus of claims 17, characterised by means (22) for trans- 
forming said reduced length filter into a minimimi phase filter. 

*25 20. The apparatus of claim 17, 18 or 19, characterised by means for 

performing the convolution step in the time domain using the discrete-time 
representation of said reduced length filter. 

21. The apparatus of claim 17, 18 or 19, characterised by means (24-36) 
30 for performing the convolution step in the frequency domain by using the 

overlapnadd method. 
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22 -fte me&od of cWm 17, 18 or 19, ch«ct«l-.d by mean, for per- 
fonning U>= convolution step in tt.e te<r«ncy domain by u»ng A. overlap- 



save method. 
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U.S. CLAIMS 

1 . A method of designing a digital filter, including the steps of 
determining a real-valued discrete-frequency representation of a desired 

5 digital filter; 

transforming said discrete-frequency representation into a correspond- 
ing discrete-time representation; 

circularly shifting said discrete-time representation; and 
applying a window to said discrete-time representation to produce a 
10 zero-padded redu<«d ler^th filter. 

2. The method of claim 1. including the step of circularly shifting said 
reduced length filter to remove leading zeroes. 
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3. The method of claim 1 or 2, wherein said real-valued discrete-frequency 
representation is formed by a noise suppressing spectral subtraction algo- 
rithm. 

4. The method of claim 1 or 2, wherein said real-valued discrete-frequency 
20 representation is formed by a frequency selective non-linear algorithm for echo 

cancellation. 

5. The method of claim 1 , wherein said window is a Kaiser window. 
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6. The method of claim 1, including the step of tiransforming said reduced 
length filter into a minimum phase filter. 

7. A digital convolution method, including the steps of 

determining a real-valued discrete-frequency representation of a desired 

30 digital filter; 

[ transforming said discrete-firequency representation into a correspond- 

ing discrete-time representation; 

circularly shifting said discrete-time representation; 
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applying a window to said discrete-time representation to produce a 
zero-padded reduced length filter; and 

convolving an input signal with said zero-padded reduced length filter. 

8. The method of claim 7, including the step of circularly shifting said 
reduced length fUter to remove leading zeroes. 

9. The method of claims 7, including the step of transforming said reduced 
length filter into a minimiim phase filter. 

10. The method of claim 7, 8 or 9, including the step of performing the 
convolution step in the time domain usix^ the discrete-time representation of 
said reduced length filter. 

11. The method of claim 7, 8 or 9, including the step of performing the 
convolution step in the fi-equency domain by using the overlap-add method. 

12. The method of claim 7, 8 or 9, including the step of performing the 
convolution step in the fipequency domain by using the overlap-save method. 

13. A digital filter design apparatus, including 

means for determining a real-valued discrete-finequency representation 
of a desired digital filter; 

means for transforming said discrete-frequency representation into a 
corresponding discrete-time representation; 

means for circularly shifting said discrete-time representation; and 

means for applying a window to said discrete-time representation to 
produce a zero-padded reduced length filter. 

14. The apparatus of claim 13, including means for circiolarfy shifting said 
reduced length filter to remove leading zeroes. 
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15. The apparatus of claim 13 or 14. wherein said window apptying means 
implements a Kaiser window. 

16. The apparatus of claim 13, including means for transforming said 
reduced length filter into a minimum phase filter. 

17. A d^tal convolution apparatus, including 

means for detennining a real-valued discrete-frequency representation 

of a desired distal filter; 

means for transforming said discrete-frequency representation into a 
correspondii^ discrete-time representation; 

means for circular^ shifting said discrete-time representation; 

means for applying a window to said discrete-time representation to 
produce a zero-padded reduced length filter; and 

means for convolving an input signal with said zero-padded reduced 

length filter. 

18. The apparatus of claim 17, including means for circularly shifting said 
reduced length filter to remove leading zeroes. 

19. The apparatus of claims 17, including means for transforming said 
reduced length filter into a minimvmi phase filter. 

20. The apparatus of claim 17, 18 or 19, including means for performing 
the convolution step in the time domain using the discrete-time representation 
of said reduced length filter. 

21. The apparatus of claim 17, 18 or 19, inchiding means for performing 
the convolution step in the frequency domain by using the overlap-add 
method. 

22. The method of claim 17. 18 or 19, including means for performing the 
convolution step in the frequency domain by using the overlap-save method. 
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ABSTRACT 

A digital convolution apparatus, includes means (12) for determining a real- 
valued discrete-frequency reptt^sentation of a desired digital filter and means 
(14) for transforming said discrete-frequency representation into a corre- 
sponding discrete-time representation. Means (16) circularly shift the discrete- 
time representation and means (18) appfy a i«dndow to the discrete-time 
representation to produce a zero-padded reduced length filter. Thereafter 
means (24-36) convolve an input signal xW with the reduced length filter. 
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(Fig. 25) 



'99 09/07 14:4 1 FAX +46 18 153050 AROS P/ArENT AB 

+46 18 

Met Patent- och reg.vert(et 
1999 -09- 0 7 

Huvudfaxen Kassan 



PATENTVERKET 3028 




Fig. 1 



Fig. 2 



Fig. 3 



li[l>-(I^l)l SCALE BY 4H] 



T T yt"! 




Fig. 4 



L+M-2 



Fig. 5 




99 09/07 14:42 FAX +46 18 153050 ARCS P/iFENT AB 

+46 18 153^8 

fcfc I Patent- ochreq-verkel 
1999 -U9- 0 7 

Hovudfaxen Kassan 2/5 



-* PATENTVERKET S|029 



1 L 



Fig. 8 



3L 



1 



L-l i 



L+M-2 



JIM 



2L i 



L+(L-1) 



L+(L+M-2) 
►^^^ 



Fig. 9 



Fig. 10 



2L *| 



2L-KL-1) 



Fig. 1 1 

il l \ " 

2L+(L+M-2) 



T 


lit T, 






, rT i 


fl .rTtr 


. , > -4 A 

OVERLAP 


^ i 


-U 





I I 



-99 09/07 14:42 FAX 18 1530 50 AROS P/kVEST AB 

+46 18 153050 



- PATBNTVERKET 12030 



Wet Patent- ochreQ-verkel 
1999 -09- 0 7 

Fig. 15 



Ilulull 



3/5 

IFFT 



I N-1 



16 



CIRCULAR 
SHIFT 



h[(ii+N/2)niodN] 



r^tTrmllTTt 



Fig. 17 



►A n 



win] 



N-1 



Fig. 18 



, K APPLY 

^ WINDOW 



MTAPS 



N-1 



«tn3=ht(tt+-N/2) modN] w[ii] 



4 


1 




O 

I 



Fig. 19 



N-1 



Fig. 21 



MAGNITUDE 



CIRCULAR 
SHIFT 

f[nM(tt*-D)inodN] 



N-1 



FFT 



1 



Fig. 20 



Fig. 22 



f[n]MIN. PHASE 



Fig. 23 



99 09/07 14:42 FAX +46 18 153050 AROS P/4rBNT AB 

+46 18 153^0 

fcfct Patent- och reg-veiket 
1999 -Ua- 0 7 



PATENTVERKET 



8)031 



Huvudfaxen Kassan 



Fig. 24 



4/5 



GET NEXT INPUT SIGNAL BLOCK 



I 



SI 



DETERMINE H[k] 



^S2 



TEIANSFORMTO 
TIME DOMAIN 



I 



CIRCULARLY SHIFT 
FILTER TO PRODUCE 
A LINEAR FILTER 



I 



APPLY WINDOW TO 
REDUCE FILTER LENGTH 



I 



CIRCULARLY SHIFT 
FILTER TO REMOVE 
LEADING ZEROES 



I 



TRANSFORM TO 
MINIMUM PHASE FILTER 



I 



TRANSFORM ZERO-PADDED 
FILTER AND INPUT SIGNAL BLOCK 
TO FREQUENCY DOMAIN 



PERFORM CONVOLUTION 
WITH INPUT SIGNAL BLOCK 
IN FREQUENCY DOMAIN 



I 



TRANSFORM CONVOLVED BLOCK 
BACK TO TIME DOMAIN 



I 



COMBINE WITH PREVIOUS 
CONVOLVED BLOCK 

1 



S3 



S4 



i^S5 



irS6 



S7 



ltSS 



irSlO 



irSll 



09/07 14:42 FAX +46 18 1530S0 AROS P/iCENT AB -> PATENTVERKET [2|032 

+46 18 153050 

Hit Patent- ochrsg-veilcBl 
1999 -09- 0 7 

Hifvudfoxen Kassan 5/5 




n 



